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Searching for Direct Stop Production in Hadronic Top Data at the LHC 
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We argue that evidence can be uncovered for stops between 300-600 GeV in 5 fb~ of 7 TeV 
proton-proton collisions if they decay into top quarks and light neutral particles. We also show that 
with 20 fb^ 1 of 8 TeV running, discovery or exclusion can be made in large regions of parameter 
space. Events with a fully hadronic top/anti-top pair and a pair of invisible decay products can be 
identified with the top-tagging of a fat jet, a single fo-tag, a missing transverse energy cut, as well 
as other kinematic cuts to reduce backgrounds with real top quarks in them. Such cuts obliterate 
the background suggesting discovery can be made with a handful of events. 



I. INTRODUCTION 

One of the central questions in particle physics since 
the 1970s has been whether or not he electroweak scale is 
"natural" [H, [1] . A theory that is not fine-tuned should 
have physics beyond the standard model that cuts off 
quadratically divergent loop contributions to the Higgs 
mass, and lives at roughly the same scale. Estimating 
using a naive energy-momentum cut off, the top one-loop 
contribution is the largest, and thus especially requires 
new light states (or other physics) to cancel its effect. In 
supersymmetric theories, it is the superpartner of the top 
quark, the scalar stop, that cancels the top contribution 
to the electroweak scale at one loop. Thus, if supersym- 
metry is responsible for naturalness, the stop should be 
accessible at the LHC, and there has been much theoret- 
ical interest in determining the near-term prospects for 
stop discovery 

Here we investigate methods of discovering stops which 
are pair-produced and then promptly decay to tops and 
long-lived neutral particles. It is natural for the stop to 
be the lightest squark if the running masses of all of the 
squarks are degenerate at a high scale. This is because 
the large top Yukawa coupling causes the stop masses 
to run down at lower energies. In addition, large mix- 
ing between the two complex scalar stops (again due to 
the large Yukawa coupling) suppresses one of the masses. 
In the same vein, the neutralino is naturally the light- 
est fermion - either the bino (superpartner of the hy- 
percharge gauge boson), or perhaps the singlino0 The 
simplest way to explain a Higgs mass between 115-125 
GeV in supersymmetry is to extend the minimal model 
to include a gauge-singlet chiral superfield, the singlino 
being the fcrmionic component. Thus, in this article, 
we explore the LHC's ability to discover a stop ('right- 
handed,' for simplicity) which decays into a neutralino, 
while all other superpartners are decoupled. 

While the LHC has not published any direct stop 
searches, the final state of tops with fix has been 



1 The neutral particle could also be a gravitino. 



searched for with 1 fb _1 of data @. This search places 
limits on the existence of fermionic top partners, but it is 
insensitive to stops because their cross section is smaller. 
This search, as well many others that would be sensitive 
to stops [7Hirjj require at least one lepton in the final 
state. Other more general searches, including, for exam- 
ple, those for fr's, jets, and $t IHj], an d those for 
jets and $r [HI El will also be sensitive to direct stop 
production, but they will certainly not be optimized. 

Below, we present a method to find purely hadronic 
stop decays with large $t using fat jets and modern 
top-tagging algorithms. Not requiring a lepton makes the 
search sensitive to the final state with the largest branch- 
ing ratio. By exploiting the properties of top decays, this 
search will be more sensitive to stops than naive searches 
with jets. Theoretical work looking at this final state 
both without [H^Tl] and with leptons fl9ll2l) studied 
the 14 TeV LHC with a much larger data set. There is 
also an analysis focusing on a lower energy LHC which 
uses leptonic channels [22j, and here we present a strat- 
egy optimized for near-term discovery with fully hadronic 
decays. 

While completing this work, a similar analysis [3| was 
published which uses many of the same methods. That 
analysis uses quantitatively different cuts than we do 
here, and it focuses mainly on 2012 LHC data, while 
we will look at both 2011 and 2012 data. 



II. RELEVANT PROCESSES 



The signal we consider is 

pp -+ it* (t X °) (f~X°) (bjjx°) (bjJX°), (1) 

where \° is & neutral long lived particle and j is a light 
flavor jet. We consider only one stop, with the second 
one potentially much heavier. We simulate the signal us- 
ing MadGraph 5 [23| and hadronize with Pythia 6 [24| . 
We use the tree level cross section from Madgraph for the 
signal and background. We do not use k-factors to esti- 
mate the next-to-leading order effects - the backgrounds 
these searches are most sensitive to are in special corners 
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of phase space, and thus k-factors may represent a poor 
approximation. 

There are two classes of backgrounds; those that in- 
clude top quarks and those that do not. Backgrounds 
which do not contain tops can be much more easily elim- 
inated with loose "top-tags," and we will discuss them 
below. The dominant backgrounds, are those with on- 
shell tops. The largest of these is tt + jets, where one of 
the tops decays leptonically producing a hard neutrino. 
A lepton veto can eliminate most of the events where 
the W decays to an electron or muon, but a r, espe- 
cially one that decays hadronically, cannot be eliminated 
in this way. We simulated background samples of ti and 
tt + 1 and 2 jets again using MadGraph 5 and Pythia 6. 
The samples were matched using the MLM procedure [25[ 
with a matching scale of 30 GeV. 

Single top + W + jets is also a significant background, 
where one of the two Ws decays leptonically. Extra ra- 
diation in these events can make it look like a tt + 
event, and we simulate t + W + 1 and 2 jets using the 
same techniques as the tt + jets sample. 

The largest background without tops is Z + jets, where 
the Z decays invisibly. Since our signal contains fully 
hadronic tops, the jets produced in association with the Z 
can in principle look like tops. As we will see in the next 
section, this does not happen very often, namely top-tags 
efficiently suppress this background. Furthermore, typi- 
cal Z's are produced with a transverse momentum ~ Mz, 
so our large cut further reduces this background. An 
even smaller background is W + jets where the W decays 
leptonically but the charged lepton is either lost or a r. 
This is small for the same reasons as Z + jets, and fur- 
ther reduced because in the W decay, the neutrino does 
not carry all of the W momentum, so an even smaller 
fraction of the events pass the J$t cut. Because we will 
require a 6-tag, we separately simulate V + bb + 1 and 2 
jets, and V+ 3 and 4 jets, where V = W,Z. 

Additional backgrounds which can produce this signal 
topology include tiZ, and diboson production, but the 
cross sections for these processes are very small and they 
need not be considered with the amount of luminosity we 
study here. 

Finally, we note that pure QCD multijet production 
can produce I^t from various sources including lcptonic 
heavy flavor decays and detector effects. The contribu- 
tion of these events can be estimated in data (see for 
example [3, HH| ) and are small for large J$t and events 
with at least one 6-tag. Furthermore, we implement a cut 
on $r/y/J2 e t as per ATLAS [13] and a cut on A0 min 
as used, for example, in this CMS study [28|. These cuts 
have negligible effects on the signal and are described in 
Section Unl 

The cross sections for all the signals and backgrounds 
at both 7 and 8 TeV energies are shown in Table Q] 



III. OUR METHOD OF SIGNAL EXTRACTION 

Our event selection is as follows. We cluster all 
hadronic activity in the event with \<q\ < 2.5 into 
"fat" jets using the Fast Jet 12911 implementation of the 
Cambridge- Aachen algorithm [3CJ,[31| with R=1.2. Then 
we take the same event and cluster "skinny" jets to mimic 
the usual experimental algorithms using the anti-fct algo- 
rithm [32j with R=0.5. We then make the following pre- 
selection cuts on our simulated samples which will have 
limited effect on the sig nal0 

• We veto all events with an isolated lepton. The 
isolation criteria is that the lepton's energy com- 
prise at least 80% of all energy in a cone of R = 0.4 
centered on the lepton. We require |ry| < 2.5 and a 
minimum px of 4 and 8 GeV for muons and elec- 
trons, respectively. 

• We veto all events that contain a hadronic "t" . We 
parameterize experimental r-tagging efficiency by 
adopting a 50% tagging efficiency for hadronic tau 
decays in which the hadronic decay products have 
p!£ d > 20 GeV and a fake rate of 2% for other jets 
|33l. l3~ij . Our tau identification uses skinny jets in 
order to match the experimental measurement. We 
assume that hadronic tau decays with less than 20 
GeV of visible pt cannot be tagged. 

• We define $t = — X); P?< where i runs over all 
visible particles and Ht — YljPTj where j runs 
over all skinny jets with pr > 25 GeV. We t hen 
require all events to have I$t I \/Ht > 5V GeV. 

• We require that the fir vector be separated by 
more than Acf> of 0.4 from each of the three hardest 



TABLE I. Tree level cross sections of signal and relevant back- 
ground processes. The number in parentheses for the signal 
rows is the stop mass. The cross section does not depend on 
the neutralino mass. In tt and single top we require exactly 
one W to decay leptonically. Here V — Z,W . 



Process 


Generator cuts 


a (fb) 


ff (fb) 




and parameters 


7 TeV 


8 TeV 


it* (340 GeV) 




254 


1.04 x 10 3 


tt* (440 GeV) 


it* -> bb + Aj + 2x 


48.8 


205 


tt* (540 GeV) 




11.8 


51.1 


tt + jets 


W-n„, p Tv > 80GeV 


16.3 x 10 3 


26.7 x 10 3 


sing, top + jets 


p Tv > 100 GeV 


4.65 x 10 3 


8.27 x 10 3 


V + bb + jets 


Z -> vv, W Iv 


1.08 x 10 3 


1.53 x 10 3 


V+ jets 


Ep t „ >80 GeV 


66.6 x 10 3 


96.3 x 10 3 



2 These pre-cuts will eliminate some signal events, but the vast 
majority of those events would not have passed our subsequent 
selection cuts. 
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leading skinny jets. Namely, min[A</>( ^T,ji)] > 
0.4. This cut quantifies the notion that I$t arising 
in QCD events tends to lie along a jet direction. 

• We require Ht > 275 GeV. This has practically no 
effect on the signal since we expect two hadronic 
top decays. 

The first two cuts reduce leptonic contributions from W 
decays. The signal, of course, also can have leptonic W 
decays and this has been used to look for stops by the 
experiments Here we are doing the orthogonal 

search of looking for all hadronic final states so we veto 
leptons. The last two cuts are designed to eliminate con- 
tributions from pure QCD, and these are cuts that the 
experiments use for this purpose. In particular, the res- 
olution on Bt in multi-jet events has been shown to be 
proportional to y/Hr [23], so requiring large Bt / \J Ht 
reduces QCD events where the Bt is mis-measured. 

The intuition for the fourth cut is that large Bt in 
multi-jet events often comes from large mis-measurement 
of the energy, but not angle, of a single jet. In that case, 
the Bt should point in nearly the same <fi direction as 
that jet, so requiring Bt to be well separated from all the 
leading jets reduces the QCD background. This variable 
is used in many hadronic searches, see for example [28j . 

After the preselection, we make our main cuts to dis- 
tinguish signal from background. 

• We require I$t > 175 GeV. 

• Among the two highest pt fat jets, we require one 
to pass a HEPTopTagger [18| . The definition of the 
algorithm is given in Appendix [XJ 

• We require the fat jet that is not top-tagged to be b- 
tagged. We assume a 70% 6-tagging efficiency and 
a 1% fake rate [H, [36|. As discussed below, we fil- 
ter the non top-tagged fat jet, and it is the filtered 
subjets that we use for 6-tagging. This procedure 
facilities a more direct comparison with experimen- 
tal procedures; it would be interesting to have di- 
rect measurements of ^-tagging in fat jets without 
using subjets. 

The I$t is the hallmark of our signal and is used to dis- 
tinguish it from events with just jets. This cut ensures 
that at least one hard neutrino is present in the back- 
ground sample (up to drastic mis-measurement, which is 
accounted for in the preselection cuts). The HEPTopTag 
helps eliminate V + bb + jets, while the &-tag reduces the 
V+ jets background. We find that requiring the 6-tag 
to be in the jet that is not top-tagged improves signal 
to background relative to allowing either fat jet to be b- 
tagged. This is because in the dominant background, the 
fat jet which is not top-tagged is often made up of exter- 
nal radiation rather than real top decay products, so this 
method of ^-tagging does a better job of suppressing this 
background. An interesting alternative to &-tagging may 




tt* 

tt + jets 

Single top + jets 
V+jets 
V+bb+jets 
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FIG. 1. Distribution of the non-HEP-tagged filtered jet mass 
for (m t -,m x ) = (340 GeV,0 GeV) at yfl = 7 TeV and C = 
5 fb _1 . This includes the preselection cuts as well as the main 
cuts, but not the transverse mass cuts. 



be to use a TV-tag on a non-HEP tagged jet, but we did 
not explore this possibility it in detail. 

After making these cuts, the dominant background re- 
maining comes from tops, specifically 



tt + jets -> (bjj) (br h v T ) + jets, 



(2) 



where is a hadronically decaying r. In order to pass 
the large Bt cut, the v T typically carries most of the 
momentum of its mother IV, although the neutrino pro- 
duced in the Th decay also has some of the momentum. 
After making these cuts, we can look at the filtered (37j 
invariant mass of the fat jet which is not top-tagged. We 
have chosen the filtering clustering radius and number of 
subjets to be Rfat — 0.3 and ns\% — 30 This variable 
is plotted in Fig. [TJ In that figure we see that there are 
some background events that pass the above cuts, but 
the mass distribution is smoothly falling. If there are 
relatively light stops and very light neutralinos, there is 
a distinct signal peak above the background around the 
top mass. Figure [5] shows that increasing the neutralino 
mass with fixed stop mass reduces the size of the signal 
peak because there is both less Bt, and the tops are less 
boosted on average, so fewer signal events pass our cuts. 
Even in Fig. [51 however, there is still modest evidence for 
signal events above the background. 

In order to further reduce background, we use two kine- 
matic variables. First, the usual transverse mass, rriT, 
which is a function of two objects, one of which could be 
the J$t of the event: 



i T (p T ,m a ,q T ,mj3) 



-2(E%E T -p T -q T ) (3) 



3 We tested pruning |3SH and trimming [39| as well, and found 
filtering to give a somewhat more peaked distribution about the 
top mass. 
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FIG. 2. Same as Fig. [T] except with (mf,m x ) — 
(340 GeV, 100 GeV). 

with 

m - \/< + (p?) 2 . (4) 

We also use the rriT 2 El[ variable which is a function 
of two visible objects as well as the J$t of the event and 
a candidate mass for the i^r, m x : 

TO| 2 (p£,m a ,p^,rna, fa T , m x ) = min 

max{mr(pr,m a , At ,%)> ™t(Pt> m /3> At i m x)l ■ 

(5) 

We have found that the optimal cuts to reduce the dom- 
inant background @ are: 

• Construct to^ 2 with the two leading fat jets and 
the missing energy of the event. Take the masses 
of the fat jets to be the measured masses and the 
missing energy to be massless, m x — 0. We require 
tot 2 > 200 GeV. This computation is done using 
the algorithm implemented by the Oxbridge MT2 
library pi]. 

• We require that mx between the J$t and each of the 
leading fat jets be greater than 200 GeV. Namely, 
min[mT(^T, j'1,2)] > 200 GeV. We take the Ifr to 
be massless, and the mass for the fat jets to be the 
measured masses. 

The tot 2 variable is designed for the events in which two 
identical massive particles decay into two identical invis- 
ible particles, as is the case for our signal process. In 
events where this is not the case, like the majority of the 
backgrounds, the rriT 2 value for the event tends be smaller 
than signal events. Therefore, this cut on tot 2 increases 
the signal to background ratio, and this fact has been 
exploited in many similar contexts, for example [TH [l8|. 

A second observation is that in the dominant tt back- 
ground (J2), the invariant mass of the hard neutrino with 




50 100 150 200 250 300 350 400 450 500 

Mj (GeV) 

FIG. 3. Same variable as Fig. [TJ but with the additional 
kinematic cuts described in the text also imposed, and the 
binning changed. 

one of the fat jets approaches (but shouldn't be larger 
than) mtop because the neutrino is carrying most the 
momentum of the W. Of course, we can only measure 
the transverse momentum of the neutrino, but the trans- 
verse mass of the neutrino with one fat jet has a cutoff 

at J mf op — m^y ~ 155 GeV. This cutoff is observed at 
parton level, but it is smeared out by hadronization and 
combinatoric background. We still find that the cut we 
implement above significantly improves signal to back- 
ground ratio. 

The same distribution as Fig.[T]is shown in Fig. [3] after 
implementing the above two kinematic cuts. Compar- 
ing Figs. [T] and 02 we see that the kinematic cuts reduce 
the number of both signal and background events, but 
they make the signal jump out above the background. 
Furthermore, the effect of these cuts on the signal point 
shown in Fig. [2] is shown in Fig. [4] While Fig. [2] shows 
only modest evidence of a signal, Fig. @] demonstrates 
that the kinematic cuts make the excess much more sig- 
nificant. The detailed results of all the cuts on the back- 
grounds and some signal points are shown in Table [TT] 

IV. CONCLUSION: ESTIMATED REACH 

We estimate the reach in the (mj, m x ) plane by look- 
ing at the invariant mass of the fat jet which is not top- 
tagged in the events where all cuts have been applied. 
This is what is shown in Figs. [3] and SJ In this sample, 
we take the events which lie in a top-like mass window 
of [150 GeV, 230 GeV] and do a simple counting exper- 
iment of the number of total events versus the number 
expected background events. We then compute a signif- 
icance assuming a Poisson distribution. The results for 
the 2011 data set of ^fs = 7 TeV and C = 5 fhT 1 are 
shown in Fig. 03 and we see that if stops exist in this 
scenario with a mass around 340 GeV along with a very 
light neutralino, then they can be discovered with cur- 
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TABLE II. Number of events passing a given cut with 5 fb 1 luminosity at 7 TeV. The LSP is assumed to be massless. 



Process 


Pre-cut 


tf T > 175 GeV 


1 top-tag 


6-tag 


m T2 > 200 GeV 


m T > 200 GeV 


it* (340 GeV) 


688 


327 


109 


50 


32 


26 


it* (440 GeV) 


150 


112 


42 


20 


16 


14 


it* (540 GeV) 


39 


33 


13 


7 


6 


6 


tt + jets 


12.5 x 10 3 


872 


248 


110 


28 


18 


Single top + jets 


1.56 x 10 3 


611 


145 


23 


8 


6 


V + bb+ jets 


906 


169 


< 1 


< 1 


< 1 


< 1 


V + jets 


9.01 x 10 3 


2.34 x 10 3 


166 


6 


3 


2 


Total Background 


23.9 x 10 3 


3.98 x 10 3 


559 


140 


39 


27 
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FIG. 4. Same as Fig. [3] but with (m t -,m x ) 
(340 GeV, 100 GeV). 



FIG. 6. Same distribution as Figs. [3] and [H but now 
with y/s = 8 TeV, £ = 20 fb -1 , and (m t -,m x ) = 
(440 GeV, 100 GeV) 
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FIG. 5. Significance of excess computed by counting events 
in rrij 2 £ [150 GeV, 230 GeV] and assuming a Poisson Distri- 
bution. This is for y/1 = 7 TeV and C = 5 fb -1 . 
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FIG. 7. Same as Fig.[5]but with yfs = 8 TeV and C = 20 fb -1 . 
Boxes with significance > 8 all have the same color. 



rent data. The same data set can also exclude stops up 
to about 440 GeV. 

Looking forward, we can repeat the analysis for a hypo- 
thetical 2012 data set with y/s = 8 TeV and C = 20 fhT 1 . 
In Fig. [BJ we see that even for larger stop and neutralino 
masses, this data set is enough to see a dramatic signal. 



We estimate the reach using the same procedure as for 
7 TeV without combining 7 and 8 TeV data sets, and 
our results are shown in Fig. [71 As expected, the reach 
improves significantly, and much of the parameter region 
up to 600 GeV can be covered. 

Our work does not take into account systematic errors 
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including uncertainties in the total cross section for the 
background. This can be ameliorated, however, by doing 
a more sophisticated shape analysis of the distributions 
in Figs. [THU rather than our naive counting experiment. 
In particular, in Figs. [3] and [4] we see that after all the 
cuts, the background is relatively smooth while the signal 
has a district peak. There will not be many events after 
all the cuts, but there is also shape information in the 
data sample without the kinematic cuts, Figs. Q] and [2] 
where there will be many more events. 

One possibility to increase the signal efficiency and ex- 
tend our reach in the (mj, m x ) plane is to decrease our 
relatively high fa cut. We have required such a high 
cut in order to reduce QCD background to a level well 
below signal. It may be possible, however, to model this 
background using data driven methods and then subtract 
it off. In this case, the thresholds on the kinematic cuts 
could potentially be increased to further reduce the top 
background and improve the search, especially for heav- 
ier neutralino masses. For example, if the stop is 250 
GeV and neutralino is massless, our current cuts pass 
less than one signal event with 5 fb _1 at 7 TeV, while 
reducing the fa cut to 100 GeV would allow as many 15 
as events to pass. 

Finally, we note that while we believe this procedure 
is optimized for finding direct hadronic stop production, 
existing experimental searches which use a similar phi- 
losophy may be sensitive to stop production. This CMS 
search [43] with 1 fb _1 uses modern top-tagging meth- 
ods to look for ti resonances. While the method used in 
this search is optimized for more boosted tops, it could 
still be sensitive to this signal. In particular, looking at 
their data sample in events with larger fa could reveal 
evidence for direct stop production. 

The experimental effort to find stops is of paramount 
importance to test whether supersymmetry solves the hi- 
erarchy problem. While this effort is underway, every 
attempt should be made to cover all different possible 
scenarios. Here we have shown that if the stop decays 
to an on-shell boosted top along with a neutralino, then 
hadronic decays of the top can be exploited to explore 
this regime with current data, and future data will allow 
for even broader exploration of the possibilities. 

Note Added: While completing this work we became 
aware of two additional theoretical works looking at sim- 
ilar searches at the LHC 0, [j| . 
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Alfred P. Sloan Foundation. DS is supported in part by 
the NSF under grant PHY-0910467 and gratefully ac- 
knowledges support from the Maryland Center for Fun- 
damental Physics. 



Appendix A: Description of HEPTopTagger 

When a boosted top decays, its decay products will 
often be collimated. This fact has been exploited to 
create several top-taggers (44T - |4r| that use various han- 
dles of hadronic top decays to distinguish them from jets 
which come from other sources. In stop searches with 
moderate stop masses such as those considered in this 
work, the rate for highly boosted tops in the final state 
is significantly lower than the overall signal rate. There- 
fore, a top-tagger which works when the tops have only 
O(l) boost would increase the signal efficiency. We find 
that the HEPTopTagger works best for modestly boosted 
tops, and that is why we choose it for this analysis. 

A detailed description of the HEPTopTagger is found 
in [18J . We give a schematic description of the algorithm 
here for completeness. The algorithm takes as input the 
fat jets described in Section lHTl It is applied to individual 
fat jets, and each one either passes or fails. The algorithm 
is as follows: 

1. Begin by undoing the clustering algorithm to look 
for subjets. When the mother jet becomes two 
daughters, if the daughter jet with larger mass has 
less than 80% of the mass of the mother jet, keep 
both jets. Otherwise, keep just the daughter with 
larger mass. Continue this process as long as sub- 
jets have mass greater than 30 GeV. 

2. Filter sub jets with resolution -Rater = 
min(0.3, Ai?jfc/2). Keep the five filtered sub- 
jets with the highest pr- If there are fewer than 
five, keep all of them. 

3. Select the set of three subjets whose combined mass 
is closest to mtop. Require that the combined pt 
of the three subjets is more than 200 GeV. 

4. Sort the three selected subjets by px, ji, J2, J3 and 
calculate pairwise invariant masses, m,\2, fTi23,mi3. 

5. If the masses satisfy 

TTLl 1 

0.2 < arctan— < 1.3 (Al) 

and 

R 

mm max • 

(A2) 

mi23 

where i? m i n = 85%mw/m top and R max = 
115% mvi// m top, then this is a top candidate. 

6. If the masses satisfy 

*4+(s) 2 ) 
<i-(^y (A3) 
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and 



7. If the masses satisfy 



W23 
™123 



> 0.35 



(A4) 



then this is also a top candidate. 



and 




< R max I 1 



m.23 



> 0.35 



"1123 

this is a top candidate as well. 



(A5) 



(A6) 



8. If the masses do not satisfy any of the three previ- 
ous conditions, then this jet is not a top candidate. 

A full explanation of why these criteria are chosen and 
what the efficiency and fake rate of this algorithm are 
can be found in [LSI ]. 
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